Scalable solid-state quantum computation in decoherence-free subspaces with trapped 

ions 
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We propose a decoherence-free subspaces (DFS) scheme to realize scalable quantum computation 
with trapped ions. The spin-dependent Coulomb interaction is exploited, and the universal set of 
unconventional geometric quantum gates is achieved in encoded subspaces that are immune from 
decoherence by collective dephasing. The scalability of the scheme for the ion array system is 
demonstrated, either by an adiabatic way of switching on and off the interactions, or by a fast gate 
scheme with comprehensive DFS encoding and noise decoupling techniques. 
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The practical accomplishment of quantum computa- 
tion requires accurate control of quantum coherent evo- 
lution to perform the information storage for quantum 
bits (qubits), the processing of information by quantum 
gates, and a means of final readout [lj. The quantum 
computer model based on ion trap systems, which en- 
codes and manipulates information via long-lived internal 
states of ions, was identified as one promising candidature 
proposal and has witnessed raoid development in the past 

decade QMM EM II H> H E3 El E3 In recent litera " 
tures [ll llfilll[ll2l |. it was suggested that quantum gates 
could be realized via certain spin-dependent Coulomb in- 
teractions. These theoretical scenarios relax significantly 
the physical constraints to execute gate operations, hence 
offer a robust manner to implement quantum information 
processing. 

The scalability of the model, i.e. the extension of quan- 
tum processors from two qubits to large numbers of ion 
units, however, is quite challenging due to the growing 
complexity of the ion vibrational mode spectrum. An 
alternative way to achieve the scalability, which is pur- 
sued by many current efforts, is to devise sophisticated 
microtrap architecture and design a series of reliable ion 
shuttling 0, @ ■ The scaling scenario concerning a large 
array of ion crystals without ion shuttling has also been 
suggested in recent proposals [llj, |l2j , where a prerequi- 
site that the gate interactions should be comparable with 
the local ion oscillation frequency is indicated. 

In this paper we propose a scheme to realize scalable 
ion trap quantum computation in decoherence-free sub- 
spaces (DFS) 1 131 with an extended unconventional ge- 
ometric scenario [lj], which possesses the advantages of 
the DFS and geometric strategies: the former is immune 
from decoherence induced by collective dephasing while 
the latter is thought to be insensitive to certain random 
errors in the operation process. We exploit the spin- 
dependent Coulomb interactions and construct a univer- 
sal set of unconventional geometric quantum gates in en- 
coded subspaces. The potential to scale up the ion array 



system for quantum computation without ion shuttling 
is further investigated. Two different interaction con- 
figurations, including an adiabatic way of switching on 
and off the interactions and a scenario to execute rapidly 
the interaction pulses combined with a noise cancelation 
technique, are presented. 

The system we employed consists of N trapped ions 
arrayed by a convenient structure. Two stationary inter- 
nal states of each ion, denoted as |0) and |1), are selected 
to represent physical qubits. We assume that the ions 
are separated by some suitable distance so that on one 
hand the ions could be addressed individually and on 
the other hand the mutual Coulomb interactions of ions 
could be taken into account. In the absence of exter- 
nal forces, the potential V of the system is normally ap- 
proximated by the second-order expansion (the harmonic 
approximation) for small vibrations around the equilib- 
rium configuration (qf\ • ■ • , q^)- The motional degrees 
of freedom of the ions is therefore treated collectively, 
and the Hamiltonian in normal coordinates reads 
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Note that the normal coordinates relate to the local one 
by Qk = Y) j DjkQji where D is an orthogonal matrix that 

diagonalizes the Hessian w y = ( d f.g ){q ( i\ ■ • ■ , gffi), 
and Wfc (k = 1, • • • , N) account for the characteristic fre- 
quencies of the normal modes. Suppose, in order to real- 
ize the qubit-qubit coupling, that two of the ions located 
at qi and qj are exerted by certain acceleration forces that 

are dependent of ion internal states by F m {t) = ffj.(t)(Xa 
with (i = The related new interaction term hence 
takes a form H F (t) = -Y^^ij f^h^a ■ We have 
assumed a general dependence of the forces on the qubit 
states, where aa with a — x,y,z represent the three 
Pauli operators acting on the states |0) and |1) of the 
qubit at site /z. It is readily seen that, by introducing 
the Fock operators a/- = (mujkQk + iPk)/ V^ muj k with 



[at, a k ] = 1, the Hamiltonian in the rotation picture with 
respect to H v n, takes the following form (setting H = 1) 

#(*) -Ebi^'+^tl^la^-^ + k., (2) 
fc=i 

where = D^f^t) for ;U = i,j and = 

To investigate the internal state evolution of the ions, 
we employ a gauged representation with respect to the 
unitary transformation G(t) = exp[— i J Q H(r)dr]. The 
evolution operator of the system is given by U(t) = 
G(t)U g (t), where U g (t) satisfies the covariant equation 
id t U g (t) = H 9 (t) U„(t) and the gauged Hamiltonian 
H 9 (t) is obtained as Ha 

H 9 (t) = G~ 1 HG — iG~ 1 dG/dt 

N 

= E4(*)^ V i j) + £ o(t), (3) 
fc=i 

where 

4(t) = / [9 h i{t)9W) + 9W)9 k M smutf - t)dt', 
Jo 

(4) 

and £o(t) is merely a c-number parameter. The equa- 
tion manifests explicitly the "spin-spin" couplings 
between the ion qubits i and j. Especially, as the 
external force is controlled with a particular configu- 
ration such that J Q T H{t)dt — 0, the transformation 
G{T) becomes an identity operator and the evolution 
operator of the system at time T is exactly the one 

U g (T) = e- 4 *( 7 >= )<T = ) = e -^HT)/2)[{a^+a^f-2] with 

= EaLi Jo J U^ dt - A ke y observation is that the 
generated transformation U g (T) contains only Pauli op- 
erators acting on the ion qubits, therefore the induced 
qubit operation is irrelevant with the ion motion degree 
of freedom, hence is insensitive to the vibrational temper- 
ature of ions. More intriguingly, similar to the case ad- 
dressed in Ref.Q, the present f/ ff (T)-gate may be viewed 
as an extended version of the unconventional geometric 
operation whose advantages have been demonstrated 
in literatures 0, H E H E3 . This sort of gate will be 
utilized as a basic one to construct the universal set of 
gate operations for quantum computation in DFS. 

We employ the pair-bit code by which the logical qubit 
is encoded in a subspace Cf as 

\0L)i = \Q)i 1 ®\l) ia , |U)i = |l)i 1 ®|0) iBJ (5) 

where i — l,---,N/2 indexes logical qubits and the 
schematic array of ions is shown as in Fig. 1. Such 
an encoding constitutes the well-known DFS [13J against 
the collective dephasing of the system-bath interaction 
Zi (8 B, where Zi = ai + ci 42 -* and B is an ar- 
bitrary bath operator. Let us denote ir x =(° J), Tr y — 



FIG. 1: Schematic of encoded logical qubits for scalable ion 
trap quantum computation in decoherence-free subspaces. 

( ° and 7T Z =( I "J as the three Pauli operators of the 
encoded logical qubit {\0l}, |li)}< A logical controlled- 
phase flip on two encoded qubits i and j, e^** w * , 
could be generated via aforementioned spin-dependent 
Coulomb interactions as below. We exploit the accelera- 
tion forces of the form F^{t) = fn(t)cr z ; , where the ob- 
jects could be the two ions {n = or alternatively 
the ions {// = 12,72} with the similar force configura- 
tion. According to the previous analysis, as the force 
configuration is designed so that at time T there are 

ij*(T)= [ T g*(t)e- i »* t dt = 0, k=l,---,N, (6) 
Jo 

the specified interactions of Eq. Q shall generate the 
transformation e - l *( T ) CT i l)cr i 3l) or e -«*(T)<r< 2) er< J2) ^ cor _ 
responding to the two different addressing of forces on 
ions /1 or respectively. Note that the evolution gener- 
ated by these interactions falls entirely into the encoded 
subspace Cf ® C? all along the time. Moreover owing 
to the simple fact of Zi = Zj = 0, the actions of both 
the quadratic operators a z Pz and o* Vi in this re- 
stricted DFS are equivalent to that of Therefore 
the gate operation e l ^ ( * ) * ( * ) could be exactly achieved by 
any of the two interacting processes described above. 

The commensurability relation (jBJ, which involves all 
of the ion oscillation modes, might be reached for large 
N systems only by an adiabatic manner to carry out 
the pushing forces on the ions. In detail, let us assume 
that fn(t) characterizing the configuration of accelera- 
tion forces is some smooth function of time and satis- 
fies <§C Wfc. Moreover, we take /^(f) undergoing 
from /^(O) = and after some finite value to an end 
point f(j,(T) = 0. It is then direct to see that there is 

J T /A«(*) e_iWfctd * = °> therefore the relations 77^ (T) = 
come into existence for all of the oscillation modes k. 
The phase of gate operation generated in this case has 
a simple form $(T) = -£ fe ( 2 M) 9l {t)g 3 {t)dt. Note 
that the currently proposed DFS scenario to implement 
qubit operations actually has already tackled partially 
the intrinsic obstacle associated with the adiabatic ac- 
tion, that the use of slow gates gives decoherence more 
time to exert its detrimental effects. 

We would also like to give some remarks on the com- 
bination of our DFS implementation of gate operations 
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with another scalable approach, i.e., the fast gate sce- 
nario by using noise cancelation techniques. In Ref. [ll| 
it was proposed that if the operation speed is comparable 
with the local ion oscillation frequency, the noise influ- 
ence due to the complexity of phonon modes could be 
significantly reduced by designing a multi-cycle configu- 
ration of kicking forces, hence the indicated scalability 
has been effectively demonstrated. From our formalism, 
as the commensurability relation © is spoiled, the trans- 
formation G M (T) (or Gfj,'(T) corresponding to the differ- 
ent ion addressing) contained in the evolution operator 
indicates actually a noise contribution, say, 

G M (T) = exp{i£ [rfr(T)a k + h.c]<rM}. (7) 

fe M={n,ji} 

It is seen that such an undesirable influence could be 
suppressed by using two cycles of force pulses with rever- 
sal configuration. The reason is that the induced noise 
effects for the two opposite evolution will encounteract 
each other by viewing that the coefficients r)„(T) and 
77k, (T) of the noise operators G M (T) and G^ (T) associ- 
ated with two circuits have reversed signs. Besides the di- 
rect extension utilizing the noise cancelation via reversed 
loop evolution, a slightly different approach to remove the 
noise effects could be achieved for our DFS gate scenario 
via combining the above specified two interactions with 
similar force configuration but different ion addressing on 
/i and //. Note the fact that in the encoded subspace of 
Cf ® C?, there are ai H) = -a^ and <^ l] = -o^ . 
The validity of the new approach relies also on the as- 
sumption that the ion array system possesses a periodic 
structure so that the relations g\ x j (f ) = gf 2 • (f) exist 
for the specified interactions regarding the translation 
invariance of the system. Since the requirement of the 
force configuration is relaxed, this method actually pro- 
vides an alternative way to implement the gate operation 
for large-scale ion array systems. 

To obtain fully the ability to perform quantum com- 
putation, one needs also construct the general rotation 
operations for single qubit units. Note that in the DFS 
schemes, nontrivial couplings between physical qubits are 
necessary even to build the gates for single logical qubits. 
We propose below a similar scheme to realize the two uni- 
versal non-commuting single-qubit gates {e l ^ x , e 1 ^ n *>} 
in the specified DFS by utilizing spin-dependent interac- 
tions. In detail, we make use of two different forces F^ a (t) 
and F'(t) (fx — 1,2) to derive the two gates e 4 "^ and 
e i4>Tr y ^ ky wmc h F^ a (t) — fn(t)ax for the former and 

= f'^t)^ with {ai 1] = a™,*™ = 4 2) } for 
the latter, respectively. The corresponding interactions 
of the form (J2J should generate the transformations 

U F {T) = e- i *^* ) ^\ U F ,{T) = e- i9 '^^y ) (8) 

at time T, respectively, provided that the relations © 
with /j, = 1,2 are satisfied. Note that in the restricted 



subspace {|0_l), |1l)}, actions of the quadratic operators 
<Jx^ offi and 0)^0^ are exactly the same as those of ir x 
and TT y . Therefore the transformations of (JSJ offer actu- 
ally the two gate operations e*^** and e 1 ^ * y , respectively. 

Protection of the state leakage throughout the gating 
period, however, needs to be scrutinized more carefully. 
It is recognized that the generated evolution of Up(t) 
and UF'{t) actually employs the ion levels out of the 
DFS, the predicted protection against collective dephas- 
ing might be spoiled during the gate operations. For con- 
venience, let us consider a simplified model with only one 
mode with frequency w involved, which accounts physi- 
cally for sideband addressing by laser beams to select 
out the particular phonon mode. We assume further a 
homogeneous dependence of the interactions on the ion 
internal states, that is, the parameters in that of @ sat- 
isfying <7i(t) = 92(t) — d(t)- Note that the evolution of 
the system now is governed by the overall Hamiltonian 

H tot {t)=H{t) + Zi®B, (9) 

where H(t) is assumed to be interactions associated with 
Ffj,(t) or F'(t) with specified parameters accordingly. It 
is readily seen that in the formerly described gauged rep- 
resentation with respect to G(f), one obtains 

Hl t (t) = H9(t) + Z i (t)®B, (10) 

where Z % (t) = G" 1 (t)Z 4 G(t) and H"(t) has a form of Eq. 
© with a degenerated parameter 

J(f) = 2 f g{t)g{t')smuj(t' -t)d£ . (11) 
Jo 

In detail, corresponding to the two interactions associ- 
ated with F M (f ) and F^ (f ) respectively, the form of Zi (t) 
could be obtained as 

Zi(t) = cosfi a (t)Z l + smfi a (t)(a ly +a 2y ), 

Z[{t) = cosr 1 a (t)Z l + smrf(t)(a ly +cr 2x ), (12) 

where the operator fj a (t) — rj(t)a + r)*(t)cv, and 77(f) = 
Jo g(T)e~ WT d,T. It is evident to see that, even for the 
evolution with perfect parameter controls with 77 (T) = 0, 
the occurrence of the final terms in the expression 112|) 
for Zi(t) and Z[{€) should inevitably mix the system and 
bath degrees of freedom, therefore spoil the desired gate 
operations. 

Notably, it happens that the decoherence effects in- 
duced above could be pined down effectively via a de- 
coupling process by devising a symmetrized multi-circuit 
evolution. For instance, the first-order of decoupling 
could be achieved via a two-cycle refocused performance 
of the interactions indicated by ?7(f + T 2 /2) = —77(f), 
where T 2 denotes the whole time period of the two cy- 
cles. In view of the relations J^ 2 rj n {t)dt = with n any 
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odd numbers and the resulted one J Q 2 smf/ a (t)dt = 0, 
one obtains readily, by using the Magnus expansion 

Ul ot {T 2 ) = f exp{-i / H 9 tot {t)dt} = e - i V*+ h *+-)T' } 

"' (13) 
the first order of the evolution with 

hi = ±- r Hf ot (t)dt 

1 2 Jo 

= I H 9 {t)dt+ j cosf) a (t)dtZi®B.(U) 

T-2 Jo Jo 

The last term of Eq. I|14|) actually contributes nothing in 
our DFS systems, therefore the decoherence effects have 
been effectively removed by the above first-order decou- 
pling process. Physically, for the two cycles with reversed 
interactions, the ions are pushed to the reverse direction, 
and the unwanted coupling of the qubits with the vi- 
brational degree of freedom induced by the dissipation 
has a reverse sign. Due to this sign reversal, the deco- 
herence effects from these two cycles encounteract each 
other. The above decoupling process could be achieved 
to arbitrary high orders by iteratively application of the 
multiple refocusing cycles |l7|. 

The expected extension of the above operation scheme 
with the comprehensive DFS encoding and noise decou- 
pling technique to the scalable system is argued as fol- 
lows. For the scenario with fast execution of interac- 
tion pulses, the refocusing concept for noises cancelation 
for the dissipative effects and for decoherence induced 
by phonon complexity is actually consistent. That is, 
we have shown that the noise cancelation by decoupling 
process is able to remove both the two kind of noises pro- 
vided that the time scale of the interaction pulse is fast 
enough comparable with the noise frequencies. On the 
other hand, for the adiabatic pushing scheme, the validity 
of the extension requires that f(t) <C uji and f(t) > r^J, 
where uji stands for the frequency of the longest wave- 
length phonon mode and denotes the relaxation rate 
of the internal state of physical ions. 

Before concluding, we would like to remark some fea- 
tures of the present scheme for a potential experimental 
implementation. Firstly, we point out that the initial- 
ization of the logical registers on \0l)i could be readily 
accomplished through laser light addressing and manipu- 
lating individually the ions i.e., initializing the states 
on |l)i 2 . Moreover, our scheme possesses the following 
advantages: (i) The noise effects of collective dephasing, 
which is reported as a major source of decoherence in the 
ion trap system ^^|, could be tackled in the scheme, (ii) 
Since the ion shuttling is not necessary in the scheme, 
one can design the ion array in any convenient geometry 
with periodic structure, (iii) There is no cooling require- 



ment since the scheme is insensitive to the vibrational 
temperature of the ions. As has been shown, all of the 
oscillation modes would not spoil the global operation 
generated by an adiabatic manner of the evolution. For 
the fast gate scenario, by making use of the noise cance- 
lation technique, there is actually very weak influence of 
the vibrational temperature on the gate fidelity [ll| • 

In summary, we have proposed a DFS scheme to im- 
plement scalable ion trap quantum computation with an 
extended unconventional geometric approach. Using the 
spin-dependent Coulomb interactions, we show that the 
universal set of quantum gates could be achieved in DFS, 
either via the adiabatic manner to switch on and off the 
interactions, or via the scenario to execute rapidly the 
interaction pulses combined with noise cancelation tech- 
niques. 
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